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Abstract 



In a previous study of the dynamics of crystals with substitutional point 
defects we had obtained simple and exact expressions for the positions of the 
perturbed crystal modes and the intensities in them in the case in which both 
the host-host and the host-impurity force-constants were taken to be central 
and nearest neighbor. Such expressions required knowledge of only one pure 
crystal lattice Green's function, the one at the defect site itself. In this paper 
we extend our previous study to incorporate non-central force-constants as 
well. We find that the same simple expressions which we had previously ob- 
tained in the nearest neighbor central force-constant case also hold without 
any modification at all in the case of isotropic force-constant changes in crys- 
tals whose host-host force-constants are isotropic, and in the particular mixed 
case in which the fractional changes in the central and isotropic components 
of the force-constants are equal to each other. In the most general arbitrary 
nearest neighbor force-constant case we obtain a reasonably compact exact 
expression for the dynamics of the perturbed modes which only involves a 
total of two pure crystal lattice Green's functions. 

I. INTRODUCTION 

Studies of nuclear resonant inelastic x-ray scattering made possible by the advent of 
dedicated synchrotron rings have generated renewed interest in the crystal impurity prob- 
lem (see e.g. [1-3]). In the typical experimental set-up x-rays are inelastically scattered off 
Mossbauer active nuclei embedded in host materials such as crystals. Such Mossbauer active 
nuclei act as impurities in the otherwise perfect host crystals into which they are inserted, 
and lead to a modification of the crystal dynamics of the host, with it being the response of 
the perturbed system rather than that of the host itself which is then measured in the in- 
elastic x-ray scattering process. The general theory for such impurity-induced modifications 
can for instance be found in [4] where exact, analytic expressions for the effect of a mass 
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change at the defect site are given. For the case of defect-induced force-constant changes 
equally exact, analytic expressions have been given [5,6] in the particular case where both 
the host-host and host-defect force-constants are central and nearest neighbor. Specifically, 
for such a situation it was found in the physically interesting central force-constant body- 
centered and face-centered cubic crystal host cases that for both of them the positions of 
the frequencies of the impure crystal were given as the solutions to 

l-p{u 2 )S{u 2 ) = , (1) 

where p{oj 2 ) is given by 

nfj*) = - 1 + Muj2 (-1 - A **M \ - M I | 2uj2 fi ^MMA (2) 
P{ ] M> + A XX (0,0)[ A' xx (0,0)) M* + u^\ A> xx (0,0)) ' {) 

M and M' are the host and defect atom masses, A xx (0,0) and A' xx (0, 0) are the pure and 
impure self force-constants at the defect site, and the function S(u 2 ) is given by 

% 2 = -1 - Mco 2 g = / dtu' 2 - 1 j (3) 

as integrated over the density of squared eigenfrequencies of the pure crystal [7] . Additionally, 
the amplitude of vibration of the defect as defined via the impure crystal defect displacement 
normal mode expansion u a (t) = (h/2u) 1 ^ 2 J2u Xa(0, 0J 2 )[a UJ e~'' lult + ale tut ] (a = x,y,z) was 
found to be given by 



|x 2 (o,c 2 )| = £l^(o,c 2 )| 



1 /M\ 2r 



MN \M' J 



[1 - p(uo 2 )S P (uj 2 )] 2 + [kuj 2 u{u 2 )p{uj 2 )Y 



(4) 



for in-band modes with u 2 < c<J r 2 nax (Sp(u 2 ) is the principal value of S(u 2 )), while being 
given by 



1*^)1 = ^ 



MKM'J [p 2 (u 2 L )T(uj 2 L ) + M/M' - [1 + p(u 2 L )Y 



(5) 



for localized modes which obey 1 — p{ui 2 L )S{uj 2 L ) = with frequencies u 2 L outside the band, 
with the function T{uj 2 l ) being given by 

For inelastic Mossbauer studies it is the quantity |x 2 (0, uj 2 ) \ which then gives the probability 
for nuclear resonant inelastic x-ray scattering at frequency cu, with the experimentally mea- 
surable PDOS (the so-called partial density of states) at energy E discussed for instance in 
[2] being given as D(E) = M'\x 2 {0, uj 2 )\v{uj) [8]. In the present paper we extend the analysis 
of [5,6] to incorporate non-central nearest neighbor force-constants as well, an occurrence 
which is also of experimental concern. 
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II. SETTING UP THE FORCE-CONSTANT CHANGE PROBLEM 



In the standard harmonic approximation the equations of motion for the displacements 
from equilibrium e~ lU!t u a (£) of the atoms of a pure 3N degree of freedom crystal lattice are 
given by 

E [A a p(e, £') - w 2 M(e')6 a p6(e, £')] u ^e) = o , (7) 

where £ ranges from to N — 1, a — x, y, z, M(£) = M is the mass of the atom at site £, and 
A a p(£,£') are the pure crystal force-constants. Similarly, for a system with a substitutional 
point impurity of mass M' located at the origin of coordinates and changed force- const ants 
A' a dl, £'), the displacements from equilibrium are given as e~ lbJ,t u a {£), with Eq. (7) being 
replaced by 

E [A a p(£, £') - w' 2 M5 a(3 5(£, £')] u p {£') = E V a/3 (£, £')u p {£') , (8) 

with the changes from the pure crystal case having been isolated in the perturbation 

V aP (e, £') = -w' 2 (M - M')6 a p5(£, 0)5(£', 0) + A a/3 (£, £') - A' af3 (£, £') . (9) 

A formal solution for the positions of the frequency modes which satisfy Eq. (8) can be 
obtained in terms of the pure crystal lattice Green's functions as evaluated in the lattice site 
representation. Specifically, one first introduces the dynamical matrix of the pure crystal 

^(k) = ^E^(0^)e-^ W (10) 

— * 

as expressed in terms of the phonon modes k of the translational invariant pure crystal, and 
then defines its eigenvectors and eigenvalues according to 

^D Q/3 (k)^(k)=^ 2 (k)^(k) , 

P 

£^(kK'(k) = <^' , E^'(k)^(k) = 5 a/3 , (11) 

a j 

and uses them to construct the pure crystal lattice Green's functions according to 

g ^ £i g) = J_ V [k)a P [ T (12) 

— * 

as summed over the three polarizations j = (1, 2, 3) and iV momentum vectors k of all the 
modes in the Brillouin zone. As constructed these Green's functions obey 

E 4*(o, QgM«>; *, = Mu 2 g a , a (u- o, O + ^E ^ m>) » ( 13 ) 

and thus immediately allow us to solve Eq. (8) in the form 
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u a (£)= £ g a ^u';e,e')v^(e',e")^(n 

t',i",f3 n 



(14) 



with the eigenmodes then being given as the solutions to the (3iV-dimensional) determinantal 
condition 

|l-G o y|=0 , (15) 

as written in an obvious notation. 

For an explicit determination of the defect amplitude of vibration |x 2 (0,u; 2 )| in any given 
case of interest, one needs to introduce the lattice Green's functions of the impure crystal 
(G) which are related to the above pure crystal lattice Green's functions (G ) via 

G = G + G VGo + G VGoVGo + ... = (l-G V)- 1 Go . (16) 

As already noted above, the frequencies of the perturbed modes are given as the solutions to 
|1 — GoV^I = 0, while the intensities of interest in the impure and pure modes at the defect 
site (respectively |x 2 (0, cj 2 )| and |Xp U rc(0' cj2 )l = 1/NM) are related by (see e.g. [6]) 

ImGU0,0) = yj Im(GoMO,0) = |x 2 (0, cj 2 )\nNu(uj 2 ) , (17) 

where v{oj 2 ) is the density of squared frequency states of the pure crystal as normalized to 
one. It is the quantity ImG xa; (0, 0) which determines the response of the system to an external 
probe, with the PDOS which gives the probability for nuclear resonant inelastic x-ray scat- 
tering at frequency uj being given by D(E) = M'\x 2 {0, uj 2 )\u(uj) = 2ujM'\x 2 {0,uj 2 )\u(uj 2 ) = 
( %jM' I nN)lmG xx (0 : 0). The defect amplitude |x 2 (0,w 2 )| modulates the response of the 
system and it is thus its determination which is needed for inelastic Mossbauer studies [9]. 

In an actual application of Eq. (17) to determine the needed |x 2 (0, cj 2 )|, the key step 
is in inverting the matrix 1 — GqV as needed to obtain G^O, 0) via Eq. (16). And even 
if one restricts to nearest neighbor force-constants only (which we will in fact do here and 
throughout) and to just a single substitutional defect, in the force-constant change case the 
matrix V a p(t,l'') will involve the defect and every single one of its host nearest neighbor 
atoms. For a simple cubic crystal for instance this defect-nearest neighbor complex has 21 
degrees of freedom (seven 3-dimensional vibrations due to one defect and six neighbors), 
while for the body-centered and face-centered cubic crystals the complexes respectively 
have 27 and 39 relevant degrees of freedom, to initially make the matrix V a p(£,£') and 
the non-trivial sector of the matrix 1 — G$V quite large. However, because of the high Oh 
symmetry at the defect site the V a p{£, £') matrix can be block diagonalized in the irreducible 
representations of the octahedral group, with the relevant decompositions in the simple, 
body-centered and face-centered cubic crystals being of the form (see e.g. [4]) 

r sc = Ai g + E g + Fig + F 2 g + 3 F\ u + F 2u 
Tbcc = A X g + Eg + Fig + 2F 2 g + A 2u + E u + 3Fi u + F 2u 

r fcc = A lg + A 2g + 2E g + 2Fi g + 2F 2g + A 2u + E u + AFi u + 2F 2u . (18) 

Then, since the displacement of the defect atom itself transforms as a 3-dimensional vec- 
tor, the defect displacements must be located entirely within the F lu modes, with all of 
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the other irreducible representations being built out of displacements of appropriate linear 
combinations which involve the nearest neighbors of the defect alone. Since it is the defect 
response to which external probes such as nuclear resonant inelastic x-ray scattering cou- 
ple, the sector of V a p(l,l') which is relevant for such scattering thus reduces to respective 
3-dimensional, 3-dimensional and 4-dimensional blocks in Eq. (16) each one of which is itself 
threefold degenerate. For the physically interesting body-centered cubic crystal the normal- 
ized 3-dimensional 3iq n mode basis is given (in body-centered cubic crystal site notation) 
by 



a 



Mo, 0,0) , 
1 



2V2 



u x (l, 1, 1) + u x (l, 1, 1) + u x (l, 1, 1) + ^(1, 1, 1) 



+u x (l, 1, 1) + u x (l, 1, 1) + u x (l, 1, 1) + u x (l, 1, 1) 
1 



«2 = T 



U 



,(1, 1, 1) + u z (l, 1,1)- u y (l, 1, 1) - u z (l, 1, 1) 



+ Uy(l, 1, 1) - U Z (1, 1, 1) - Uy(l, 1, 1) + U Z (1, 1, 1) 

+Uy(l, 1, 1) + U g (l, 1, 1) - Uy(l, 1, 1) - U Z {1, 1, 1) 
+Uy(l, 1, 1) - U Z (1, 1, 1) - Uy(l, 1, 1) + U Z (1, 1, 1) 



(19) 



while for the equally interesting face-centered cubic crystal the normalized 4-dimensional 
4F iM mode basis is given (in face-centered cubic crystal site notation) by 



a 
ai 



^(0,0,0) , 

1 



2V2 



u x (l, 1, 0) + u x (l, 1, 0) + ^(1, 0, 1) + u x (l, 0, 1) 



«2 



a 3 



+u x (l, 1, 0) + u x (l, 1, 0) + u x (l, 0, 1) + ^(1, 0, 1) 

1, 0) + u y (l 1, 0) + u z (l, 0, 1) + u z (l, 0, 1) 

-Uy(l, 1, 0) - Uy{l, 1,0)- U z (l, 0, 1) - U Z (1, 0, 1) 

L l u x (0, 1, 1) + u x (0, 1, 1) + u x (0, 1, 1) + u x (0, 1, 1) 



(20) 



Degenerate with each of these bases are two others, one based on 1^(0,0,0) and the other 
on u z (0,0,0). With the body-centered and face-centered cubic crystals being each other's 
reciprocal lattice, in the harmonic approximation where momentum and position are treated 
equivalently, the defect responses in the two cases must be identical [5], and it will thus suffice 
in the following to treat just one of the two of them. And with the body-centered cubic basis 
of Eq. (19) having lower dimensionality than the face-centered cubic crystal basis given in 
Eq. (20), in the following we shall treat the body-centered cubic alone. (While we do not 
treat the simple cubic crystal case here, for it one can anticipate results analogous to those 
we shall provide below for the body-centered cubic case.) 

In terms of the basis of Eq. (19) the matrix of the pure crystal Go in the body-centered 
cubic Fi u mode can be written very compactly as [6] 
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/ g 2^2 9l 4g 2 \ 
[G ] Flu = 2^2 9l Q y/2R\ , (21) 
V Ag 2 ^/2R S + T) 



where we have introduced the notation 

2o = fe(000) , ^1 = ^(111) , ^ 2 = ^(111) , R = g xy {222) + g xy (220) , 

Q = 9o + fe(222) + g xx (200) + g xx (022) + 2g xx (220) + 2^(020) , 

S = go + g y y(222)-g yy (020)-g yy (202) , T — g yz (222) — g yz (022) (22) 

which takes advantage of the translational invariance of the pure crystal lattice to set 
g a p(u; £, £') = g a p(w, £ — £', 0) = g a p{£ — £'). The pure crystal Green's functions which appear 
in Eq. (21) are not completely independent of each other, as some of them are related via 
the general Eq. (13). Specifically, if Eq. (13) is restricted to nearest neighbor force-constants, 
for the body-centered cubic we obtain the relations [6] 

A xx {0, 0)g + 8^(111)01 + 16^(111)^ = 1 + Mu 2 g , 
84^(111)02 = , 

A xx (0, 0)^1 + A XX (U1)Q + 2^(111)^ = Mu 2 gi , 

A xx (0,0)g 2 + A xx (lll)R + A xy (lll)(S + T) = Mco 2 g 2 . (23) 

(Our notation here is to set A a p(£,£') = A a p(£ - £' , 0) = A a p(£ - £'), but to use A a p(0, 0) 
to denote A a p(£ — 0, £' — 0).) In establishing Eq. (23) we have used the fact that for the 
body-centered cubic crystal the lattice delta function (1/N) J2n e ik ' R ^') with R(^') = (p, q, r) 
has the property that it is equal to one if p + q + r is an even integer and equal to zero 
otherwise [10]. As we see, via the use of Eq. (23) only two of the Green's functions which 
appear in Eq. (21) are independent. Evaluation of Eq. (16) can thus involve no more than 
two independent pure lattice Green's functions. 

In order to be able characterize the various pure crystal force-constants A a p(£,£') which 
appear in Eq. (23) as well as their impure A' a/3 (£,£') counterparts, we recall that in terms of 
the two-body host-host interatomic potential 4>(r), the force-constants between a host atom 
vibrating around site R a {£) and one vibrating around the origin are defined as 

A (£0)- ( \ - K W ^'^)r(£)R(£) 0,(r) * (24) 
Aa ^ 0) " " \du a {£)dup{£) Jo"" [— ~ ) R ^ R ^ - ' (24) 

as calculated at the equilibrium lattice separation between the host atoms. Force-constants 
for which <f>'(r) just happens to vanish when the atoms are in their equilibrium positions 
are referred to as being central, while those associated with potentials which obey <f>'(r)/r = 
0"(r) at equilibrium are referred to as being isotropic. In terms of the various A a p(£, 0) 
with £ 7^ the self-force-constant at the origin is then given via Newton's third law as the 
summation 

A a p(0,0) = -J2A a p(0,£) . (25) 

For a nearest neighbor pure crystal the force-constants can thus be characterized in terms 
of two parameters, viz. 
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(26) 



in terms of which we obtain 



f a + (3 a 



a 



a + (3 —a —a 



a a + (3 a , A^lll) = —a a + (3 a 



\ a 



a a + (3 



—a 



a a + (3 



and 



(27) 



A xx (0,0) = -8(a + (3) , A xy (0,0)=0 . 



(28) 



Given Eqs. (27) and (28) we can now rewrite the Green's function relations of Eq. (23) in 
the convenient form 

8(a + p)g x = 8(a + P)g + 1 + Mu 2 g , 
92 = , 

(a + (3)Q = Mu 2 g! + 8(a + (3)gi - 2aR , 

a(S + T) = -(a + (3)R . (29) 

When a defect is introduced substitutionally at a lattice site, the interatomic potential 
<p(r) between it and a host neighbor will in principle differ from that between pairs of host 
atoms, leading to modified force- const ants 



' d 2 0(r) \ 
du a (i)dup(£)J 



and 



^(o,o) = -E^(o, 



(30) 



(31) 



as again calculated at the host atom lattice equilibrium separation (the defect being in- 
serted substitutionally), with the modified force-constants then being characterized by two 
parameters 



a — — - 



(r)- 



<f>'(r) 



4>'(r) 



(32) 



whose relation to a and (3 can be arbitrary [11]. For the impure system we thus obtain 



4*(ui) = 



(a + (3 



a 



a 



\ 



a a + (3 a 
\ a a a + (3 J 



, 4*(111) = 



/ a + (3 —a —a \ 

—a a + (3 a 
\ —a a a + (3 ) 



(33) 



and 



A' xx (O,O) = -8(a + 0) , AL (0,0)=0 



(34) 
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Finally, on introducing the parameters 

X = a-a , Y = P-j3 . (35) 
we can write the matrix elements of the perturbation V a p(t,£') in the F lu mode as 



[V]f 1u = 



f -u 2 {M - M') - 8X - 8Y 2^2{X + Y) AX \ 

2^2(X + Y) -(X + Y) -y/2X 

V 4X -y/2X -2X-YJ 



(36) 



Armed with Eqs. (21), (29) and (36) we can now proceed to an evaluation of (1 — G V) 1 G . 



III. SOLVING THE FORCE-CONSTANT CHANGE PROBLEM 



To solve the problem we first proceed symbolically and set 



[1 " G V] Flu 



(a b c\ 

d e f 
\g h i ) 



where 



a 
b 

c 
d 

e 

f 
9 
h 

i 



l-u; 2 (M'-M)g + 8(X + Y)(g -g 1 ) , 
-2y/2(X + Y)(g - gi ) , 
-AX(g - gi ) , 

-2 v /2[cj 2 (M' - M)gi - (X + Y)(8 gi - Q) + 2XR] 
1-{X + Y)(8 gi -Q) + 2XR , 
y/2[-X(B gi -Q) + (2X + Y)R] , 
— 4(X + Y)R — AX{S + T) , 
V /2[(X + Y)R + X(S + T)] , 
1 + 2XR+ (2X + Y)(S + T) . 



In terms of these symbolic quantities the inverse matrix is given by 

ei — fh he — ib fb — ce\ 



[1 - G V]- F l 



fg — di ia — gc cd — af 



dh — ge gb — ha ae — bd J 
where the determinant of [1 — GoV]f 1u is given by 

A = a(ei - fh) + d(hc - ib) + g(bf - ce) . 



(37) 



(38) 



(39) 



(40) 



From Eqs. (16) and (21) and recalling the vanishing of g 2 , it then follows that the defect 
site component of the impure crystal Green's function G is given in closed form by 



G xx (0, 0) = -[(ei - fh)g + 2^2^ - ib) 9l ] 



(41) 
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with the matrix 1 — GqV having been inverted analytically. 

For computational purposes we note that in both A and G xx (0, 0) it is just three symbolic 
combinations, viz. ei — fh, he — ib and bf — ce, which are needed. From Eq. (38) these 
combinations are readily found to evaluate to 



ei 



fh 



Q) 



1 + AXR + (2X + Y)(S + T)-(X + Y)(8 gi 
+ (3XY + Y 2 )[-(S + T)(8 gi -Q)-2R 2 } , 
hc-ib = 2^/2( 9o - gi )[X + Y+ (3XY + Y 2 )(S + T)} , 
b f- ce = -4( go - gi )[-X + (3XY + Y 2 )R} . 

Further algebra shows that the numerator in Eq. (41) evaluates to 

NUM = (ei - fh)g + 2^j2(hc - ib) gi 

= g [l + AXR + (2X + Y)(S + T) + (X + Y)Q\ 

-S(X + Y)g\ + (3XY + Y 2 )[(S + T)( 9o Q - 8g 2 ) - 2g R 2 ] , 

while the denominator evaluates to 

A = 1 + AXR + (X + Y)(8g - 16 gi + Q) + (2X + Y)(S + T) 

+ (3XY + Y 2 )[(S + T)(8g - 16 gi + Q) - 2R 2 ] - 8(M - M')uo 2 (X + Y)g\ 
+ (M - M')uo 2 g Q [l + AXR + (2X + Y)(S + T) + (X + Y)Q\ 
+ (M - M')u 2 {3XY + Y 2 )[{S + T)(g Q - 8gf) - 2g R 2 } . 



(42) 



(43) 



(44) 



To simplify the problem further we now utilize the pure crystal Green's function rela- 
tions given in Eq. (29), to find, following a fair amount of algebra, that the numerator and 
denominator in G xx (0,0) reduce to 

9o 



NUM 



4 91 



and 



a(a + (3) 

[1 + (M - MVgp] 
[1 + Muj 2 go - M'iu 2 gi ] 



a(a + 13)- 2f3 2 XR - (3a 2 + 2a[3 + [3 2 )YR 
a + (3){3XY + Y 2 )R - a(X + Y) 



(45) 



a(a + (3) 



a(a + P)- 2j3 2 XR - (3a 2 + 2a(3 + (3 2 )YR 
(a + (3)(3XY + Y 2 )R - a(X + Y) 



where 



9i = 9o + 



[1 + Mu 2 9o } 



8(a + (3) 

In terms of the convenient functions 

(a + (3)(3XY + Y 2 )R (a + f3)((3 - (3)(3a -3a + f3-(3)R 



(46) 



(47) 



R 



a(X + Y) 
2((3X - aY) 2 



a(a — a + (3 — (3) 
2(a(3-(3a) 2 



(a + (3)*(3XY + F 2 ) ( a + /5) 2 (/3 - (3)(3a - 3a + (3 - (3) 



(48) 



9 



which obey 



2(a(3-(3a) 2 R 



a(a + P)(a — a + (3 — (3) 
the numerator and denominator in G xx (0,0) can be written more compactly as 



(49) 



NUM 



A xx (0,0) 



go + 



(1 + Mu 2 g ) ( A xx (0,0) 
A xx (0,0) U^(0,0) 



- 1 



Mu 2 



A xx (0,0) 



M 



p(<S)S{tS) - —S(u 2 ) - 1 



1 -R 
1-R 



+ p[S(u 2 ) 



l\R 



(50) 



and 



+ 



M'u 



I, ,2 



M'uj 



I, ,2 



MM'u^go MM'u 4 g 



A xx (0,0) A' xx (0,0) 



R 



A xx (0,0) A' xx (0,0) 

— jiR [l + (M — M')uj 2 g Q 



M' j A' xx (0,0) 
M 1^(0,0) 



P (lu 2 )S(uj 2 



R 



fiR 



(51) 



where piyj 2 ) and S(u 2 ) = —1 — Mu 2 go were introduced earlier in Eqs. (2) and (3). The 
form for G xx (0,0) = NUM/A implied by Eqs. (50) and (51) is our main result, an exact 
relation for a body-centered cubic crystal under the sole assumption of nearest neighbor 
force-constants, with there being no restriction on the strengths of the central and non- 
central force-constants in the pure crystal or the amount by which they might change in the 
presence of the defect, with the evaluation of G xx (0,0) requiring a knowledge of only two 
pure crystal lattice Green's functions, go and R. 

While the form of Eqs. (50) and (51) is completely general, great simplification occurs 
whenever we can set pR = 0, i.e. whenever the force- const ant changes obey 



since then G xx (0,0) reduces to 
GW0,0) = 



af3 — j3a = , 



[ P {lo 2 )S{u 2 ) - {M/M')S{u 2 ) - 1] 
M'uj 2 [1 - p(u 2 )S(uj 2 )} 



(52) 



(53) 



to then only depend on one pure lattice Green's function alone, the one at the defect site 
itself. In this restricted case we see that the positions of the eigenmodes are given as the 
solutions to 



l-p^ 2 )^ 2 ) = , 
while on giving u 2 a small imaginary part in the complex plane and recalling that 



1 

uj' 2 + it) 



L_^-i^ 2 -J 2 ) , 



(54) 



(55) 
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find that the imaginary part of G xx (0, 0) evaluates to 



[1 - P {u 2 )S p {uj 2 )] 2 + {kuo 2 v{uo 2 ) P {uj 2 )\' 



(56) 



With the imaginary part of the pure crystal go evaluating to Img = 'ku{uj 2 )/M = 
N\xl urc (0,u 2 )\nv(cu 2 ), Eqs. (1) and (4) thus follow [12]. With Eqs. (1) and (4) having pre- 
viously been obtained under the central force-constant restriction, we see now that they in 
fact have far greater validity, since as well as the pure central force-constant solution to Eq. 
(52) in which (3 — (3 — 0, Eq. (52) also admits of other solutions. Specifically it admits of a 
pure isotropic force-constant solution in which a = a = 0, and also of solution in which the 
central and isotropic components of the force- const ants undergo the same fractional change 
a/a = (3/(3. 

While the family of solutions which obey a/3— (3a = encompasses a large and interesting 
class of force-constant changes, to go beyond this set of solutions requires including the 
second Green's function R. Unlike the Green's function g whose evaluation depends only 
on the density of states of the pure crystal according to 

go = g xx (u;0,0) = ^[g xx (uj; 0, 0) + g yy (u; 0, 0) + g ze (u; 0, 0)] 

i [^'(k)^(k) + ^(k)^(k) + ^'(kH(k)] 



3iVM^ K(k)-^] 

as can be seen from the definition of Eq. (12), any other pure lattice Green's function such 
as R will involve a sum over momentum modes which cannot be reduced to an equivalent 
sum over eigenfrequencies. However, as noted in [5] it is possible to simply the evaluation of 

— * — * 

the lattice Green's functions which involve a non-zero R(€') — R(£) to some degree, as it is 
possible to perform the sum over polarizations analytically. Specifically, it was shown that 
the general g a p{tjJ\ i, 0) with £ ^ can be written in closed form as 

9afi (u; £, 0) = £ X af3 (k)e-^ (58) 

k 

where the different components of A Qj a(k) can be written in terms of the pure crystal dy- 
namical matrix D a/3 (k) of Eq. (10) as 

A f * \^-uj 2 {D yy + D zz )+D yy D zz -D 2 yz } 

Xx[ } [^-^ED xx + u 2 (ED yy D zz -ED 2 xy )-UD xx + ED xx D 2 z -2UD xy ] 1 ] 



and 



\^ D xy D xy D zz -\- D xz D yZ j 



, L ^ xy ^xy-^zz i J ~'xz j ~"yz\ (60) 

xy[ } = K - cu 4 E D xx + ^ 2 (E D yy D zz - E D 2 xy ) -UD XX + E D xx D 2 yz - 2 n D xy ] [ } 

(the sums being taken cyclically), to thus reduce a determination of the pure lattice Green's 
functions to a straightforward sum over the Brillouin zone. 
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[7] In Eq. (2) we have used the relation ou 2 mx = 2A XX (0, 0)/M which holds [5] for any pure 
harmonic cubic crystal in the nearest neighbor force-constant approximation. In Eq. (3) 
we have introduced the pure crystal lattice Green's function g = g xx {uj\ 0, 0) which we 
will describe below. The pure crystal density of squared frequencies v{uo' 2 ) = v(uj')/2uj' 
which is associated with go is normalized to J^ max du l2 v{oo 12 ) = 1. 
[8] The PDOS is referred to as being partial since not all the nuclei in a given unit cell are 

necessarily Mossbauer active. 
[9] In general for a perturbation Q the response function is given by K{uj) = ulmTiGQ 
where the trace is taken over the relevant Oh representations. 

[10] The four relations given in Eq. (23) respectively follow by setting (£', a, a') equal to 
(0,x,x), (0,x,y), (lll,x,x) and (111, rr, y) in Eq. (13), restricting the sum on t to 
the origin and its eight first neighbors, and using the symmetry relations A xy (0,0) = 
0, A xx {\\\) = A xx (lll), A xy (lll) = -4^(111), ^(0,0) = 0, fc(Hl) = fc(lH), 
<^(111) = -g xy (lll), g xy (002) = g xy (020) = g xy (200) = g xy (022) = g xy (202) = 0. In 
passing we note that in Ref. [6] the 8A xy (lll)g2 = relation (a relation which entails 
that the pure body-centered cubic lattice Green's function g 2 is actually zero in the 
nearest neighbor approximation) was not considered as it was not needed there, with 
the conclusions reached in the pure central force-constant study of [5,6] not in any way 
being dependent on the value of g-i since the use of the three other Green's functions 
relations given in Eq. (23) still allowed for its elimination in favor of the defect site go 
anyway. However, for the simplification of the non-central force-constant calculations 
which are presented in this paper, the vanishing of gi turns out to be quite crucial. 

[11] Even if the pure crystal lattice separation is such that nearest neighbor atoms just 
happen to be at the same distance as the minimum in the two body host-host interatomic 
potential (central force- const ant crystal case), it does not follow that the minimum in the 
two body host-defect interatomic potential would occur at precisely that same distance. 

[12] In passing we note that there is a difference between the usage of the real uj 2 equation 
1 — p(uo 2 )S(uo 2 ) = of Eq. (54) which gives all the modes of the impure crystal and 
the usage of the complex uo 2 condition 1 — p(u 2 )Sp(uj 2 ) = at which ImGr xx (0, 0) 
of Eq. (56) is resonant, with there being at best only a handful of frequencies for 
which this resonant condition is satisfied, and possibly none at all if p{oo 2 ) is such 
that the 1 — p(u 2 )S(u> 2 ) = equation has a localized mode solution outside the 
pure crystal frequency band. (Outside the band there is no need to distinguish be- 
tween S(u 2 ) and Sp(u 2 ) since S(u 2 ) has no poles there.) Specifically, in applying 
Eq. (54) to identify all of the allowed impure crystal eigenmodes one uses the dis- 
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crete form of S(u 2 ) given according to Eqs. (3) and (12) and Eq. (57) below as 
S^ 2 ) = _i _ Muj 2 g xx (u; 0, 0) = -(1/N) ^H(^ 2 (k)]/[cu 2 (k) - uo\ to find 

that the modes of the impure crystal lie either between the modes of the pure crystal 
where S{oj 2 ) diverges, or below the lowest pure crystal frequency mode or beyond the 
largest pure crystal frequency mode, with the discretized spectra of the pure and impure 
crystals having no common eigenfrequency. (Despite being a sum on both k and j, the 
<7* J '(k)<7-£(k) term in the expression for S(uj 2 ) only permits the j sum to contribute once, 
since in the basis in which D Q( g(k) is diagonalized we can choose the polarization vec- 
tors cr^(k) so that only one polarization eigenvector has a non-zero x component. The 
function S(u 2 ) thus only possesses N poles, with the full 3N dimensionality of the eigen- 
frequencies then being recovered via the additional eigenmode equations for 1^(0, 0, 0) 
and u z (0, 0, 0).) This lack of any common eigenfrequency between the pure and impure 
systems holds no matter what values one uses for the M/M' and A xx (0, 0)/A' xx (0, 0) 
parameters which appear in the 1 — p{uj 2 )S(uj 2 ) = condition (H. J. Lipkin and P. 
D. Mannheim, Bounds on Localized Modes in the Crystal Impurity Problem, cond- 
mat/0510542), with the threefold degenerate 1 — p(uo 2 )S(uo 2 ) = condition always gen- 
erating a total of 3N perturbed crystal eigenfrequencies. For perturbations which obey 
the af3 — j36i = condition of Eq. (52) then, the threefold degenerate F\ u sector condi- 
tion 1 — p(o; 2 )S'(a; 2 ) = generates the full set of 3iV impure crystal eigenmodes, with 
there being no possibility of any further modes being generated via eigenvalue condi- 
tions associated with any of the other irreducible Oh representations that are given in 
Eq. (18). 
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